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Abstract 

We have extended the classical over-barrier model to simulate the neutral- 
ization dynamics of highly charged ions interacting under grazing incidence 
with conducting and insulating surfaces. Our calculations are based on sim- 
ple model rates for resonant and Auger transitions. We include effects caused 
by the dielectric response of the target and, for insulators, localized surface 
charges. Characteristic deviations regarding the charge transfer processes 
from conducting and insulating targets to the ion are discussed. We find 
good agreement with previously published experimental data for the image 
energy gain of a variety of highly charged ions impinging on Au, Al, LiF and 
KI crystals. 
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I. INTRODUCTION 



Within the past decade, a rapidly increasing number of research projects has been de- 
voted to the investigation of interactions between highly charged ions (HCI) and surfaces 
(for recent reviews see These activities are of importance for present and future appli- 

cations, such as semiconductor fabrication, nanostructure technology, and surface chemistry. 
They are also of interest to basic research due to the challenging interplay of fundamental 
electronic interactions to be considered in the detailed understanding of the highly complex 
interaction dynamics. By now, a certain level of consent has emerged with regard to charge 
exchange and ionization processes that take place before an HCI gets in close contact with 
a metal target surface, and a mainly classical approach, the "classical over-barrier model" 
(COM), first presented by Burgdorfer, Lerner, and Meyer 0,^, was found to adequately rep- 
resent the most important physical aspects of the electron capture, recapture, and emission 
sequence [pHTlH. 

Typically, an incident HCI captures several conduction band electrons at large distances 
from the surface into highly excited states which leads to the temporary formation of a 
"hollow ion" in front of the surface. At ion-surface distances that are smaller than or about 
equal to the classical radii of active HCI orbitals, the theoretical description becomes more 
difficult due to the strong perturbation of the initial electron distribution of the surface and 
the intricate molecular dynamics involved. For this reason, most first-principle calculations 
have been applied to incident ions in low charge states |l^-[l5|. For higher incident charge 
states the detailed quantum mechanical treatment is complicated by a large number of ionic 
states that are energetically degenerate with the target conduction band, and a first-principle 



approach remains a formidable task |[T6|-[19|. 

Most experiments with incident HCI have been performed for conducting and semicon- 
ducting surfaces with typical workfunctions of about 5 eV. These experiments focused on 
total electron yields |2^, energy-resolved Auger pTj-p^ and X-ray p5| , |26[| spectra, as well 



as defiection angles p^, and ion-neutralization [p^-|30[] measurements. Recently, several 
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experiments have been carried out, where HCI beams are incident on insulating surfaces, 



primarily on ionic crystals, such as LiF pT|-[3^. The unique band structure of LiF with a 



large workfunction of 12 eV and a wide band gap of 14 eV (Fig. |1]) that elevates the an- 
tibinding 2p-band above the vacuum level provides an interesting opportunity to scrutinize 
previous theories about the role of the conduction band in the charge exchange process. 
In contrast to metal surfaces, LiF and other insulating surfaces do not provide unoccupied 
conduction band states into which resonant loss from excited projectile states might occur, 
and pronounced differences are expected in the neutralization dynamics of HCIs in front of 
metal and insulating surfaces. Furthermore, the capture of electrons from an insulator leads 
to the local accumulation of positive surface charges that modify the charge transfer dynam- 
ics in comparison with metals. The capture-induced accumulation of localized charges on 
insulator surfaces has very recently been addressed in a few independent theoretical stud- 
p5|^7| and is a central aspect of the extentions to the COM discussed in this paper. 



les 



Measurements on LiF surfaces [^] clearly exhibit the expected discrepancies with respect 
to conducting targets and can coherently be interpreted by a retarded onset for electron 
capture and characteristic deviations in the succeeding charge transfer sequence. 

In this paper, we adopted the basic framework of the COM suggested by Burgdorfer et 
al. for the interaction of slow HCIs with metal surfaces Different versions of COMs 

have been applied to successfully model charge exchange, energy gain, and trajectory effects 
in collisions of HCIs with atoms and clusters ( Ceo) M,!!!- Within the COM, the 



neutralization dynamics is described by means of an effective single electron potential which 
governs the classical motion of electrons that are going to be either resonantly captured into 
hydrogenic projectile levels or resonantly lost to unoccupied states of a metal target con- 
duction band. For charge transfer to occur, an active electron must overcome the potential 
barrier between projectile nucleus and target surface. Electron transfer becomes classically 
possible if the electron initially occupies a state that lies above this potential barrier and if 
vacancies exist in the resonant final state into which the electron transits. As the projectile 
moves along a classical trajectory, both projectile and target levels experience variable level 



shifts and change their relative energetic positions with respect to the potential barrier, the 
position and height of which also changes as a function of the projectile-surface distance. 

In order to perform simulations involving insulating surfaces we have extended the orig- 
inal COM by modifying the dielectric response of the surface to the external projectile 
charge and by including local surface charges. These local charges are built up on the in- 
sulator surface during the charge transfer sequence and decay on a time scale that is given 
by the conductivity of the insulator. The local surface charges influence active electrons 
and the projectile motion. In this work we will discuss the influence of these excess surface 
charges on electron transfer and projectile deflection in detail and comment on the recent 
and related theoretical work of Borisov et al. and Hagg et al. |]36 |. 

We have organized this paper as follows. In Section |I| we review the main physical 
elements of our COM simulations, such as effective potentials (Sec. |II A|) , local workfunc- 
tion changes (Sec. [II B|) , electronic transition rates (Sec. [II C|) , and the projectile motion 
(Sec. [II D[ ). In Section [III A[ we discuss in detail the neutralization dynamics in front of 
the surface in terms of the evolution of level occupations, potentials, projectile charge and 
motion above metal and insulator surfaces. In Section [llllj[ we compare our results with 
previously published experimental and computed data on image energy gains of the projec- 



tile over a wide range of initial charge states. Our conclusions are contained in Section |V 
We use atomic units (a.u.) throughout this paper unless specified otherwise. 



II. OUTLINE OF THE EXTENDED COM 

A. Potentials seen by an active electron 

Within the dynamical COM charge exchange is described in terms of classical charge 
currents between energetically shifted valence states of the target and shifted hydrogenic 
projectile levels. These continuous charge currents correspond to electronic transition rates 
for resonant capture and loss and occur as soon as the potential barrier Vb of the total 
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effective potential Vtot drops below the target workfunction W. The total potential acting 
on an active electron is given by 

Vtot{<lp, X, Z, Xp, R, t) = Vproj{(lp, X — Xp, Z — R) + Vim,p{qp, X — Xp, z + R) 

+ Viyn,e{z) + Vlocal{x, Z, t) (2.1) 

where x and z will denote the electronic coordinates in the collision plane parallel and 
perpendicular to the surface-projected motion, respectively. The projectile distance from 
the surface is denoted by R. The projectile coordinate along the projection of the trajectory 
on the surface is Xp. The coordinate xq < refers to the location on the surface where, 
at time to? charge transfer starts (Fig. |]). The origin of our coordinate system is located 
on the intersection of the topmost lattice plane (at z = = 0) and the collision plane. 
The jellium edge is located half a lattice constant above the uppermost lattice plane of the 
crystal. The potential saddle is located at Zh- The projectile is assumed to reach its point of 
closest approach to the surface at time t = 0, and the surface projection of this point defines 
X = Xp = 0. For our applications in this paper, it is sufficient to consider trajectories with 
surface projections along the [100] direction that intersect surface lattice points and define 
the projectile coordinate Yp = (see Section IIIB2 below). A more general approach would 
average over many trajectories with Yp coordinates inside a surface unit cell. 

The first term in (|2.1|) represents the interaction of the active electron with the projectile 



and is modeled by the Coulomb potential 



VproAlp, x-Xp,z-R) = - (2.2) 

[x - XpY + [z - Ry 



The projectile charge 

qp{R) = qnuc,p - Yl (^n{R) (2.3) 

n 

depends on the nuclear charge qnuc,p of the HCI and the projectile shell occupations a„. The 
index n labels the principal quantum number of projectile shells. 

Surface charge distributions produced in response to the external charges of the projectile 
and the active electron are included in (|2.1|) in form of the projectile image potential Vim,p and 



the self-image potential of the active electron Vim,e- These potentials can be derived within 
linear response theory. Along an axis that is perpendicular to the surface and includes the 
projectile nucleus (i.e. for x = Xp), approximate expressions for these potentials, appropriate 
for grazing-incidence collisions, are given by |^ 



^im,p(gp, 0, 2: + i?) = / duRe\ .Ko\-{z + R-2zi^)\ (2.4) 

and the self-image potential of the active electron 

VimA^) = -— TdujRe f ^^441 • ^0 f2-(^ - ^^rn)] (2.5) 

nvp Jo \l + e{LU) J \ Vp J 

where Vp denotes the projectile velocity and Kq is a modified Bessel function. 

Equations (p.4|) and (|2.5|) refer to the dielectric response of the target material to a moving 



external charge in the undispersive approximation [40|, for which the dielectric function 



e{k,uj) is independent of the momentum k. Following reference we approximate the 
dielectric function 

e{u) = eoo + . / X (2-6) 

by its static (eo = e(0)) and optical (eoo = e(oo)) limits, the resonance frequency uq, and a 
positive infinitesimal constant 7. Table I contains these constants for the two ionic crystals 
used in this work. Both image potentials are referred to the image plane at > 0. In 
applications to metals, we identify the image plane with the jellium edge, such that Zim 
becomes equal to half a lattice constant. For ionic insulator crystal targets, we inserted 
the negative ion radius for Zim, i.e. we assume that the induced positive image charge is 
located in the vicinity of the high density anionic electron cloud above the target. A similar 
independent study of the dynamic dielectric response of alkali halides was recently published 
by Hagg et al. ||36| in which, as in (p.4|), (p. 51), and (|2.6|), the linear response theory result of 
Abajo and Echenique was combined with a non-dispersive single-pole fit to the dielectric 
function. The main difference to our approach appears to be the use of a small but finite 
damping constant 7 in the work of Hagg et al.. We also note that Barany and Setterlind 
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0] have developed a COM for capture from a dielectric sphere of radius a and frequency 
independent, non-dispersive dielectric constant e. Their limit a ^ oo corresponds to an 
insulating surface with simplified image interactions that include dielectric screening effects 
in terms of a frequency independent, multiplicative factor (1 — e)/(l + e). 

For metal surfaces, we can take the limits e ^— > oo in (^.41) and (|2.5|) and obtain the simple 



asymptotic forms 



VX^'iqp, 0,z + R)= (2.7) 

1 



F.-'^f'(^) = ^. (2.8) 

In order to avoid the unphysical singularity of the electron self-image potential at z = Zim, 
we truncate and steadily connect the total image potential to the bulk potential given by 
the lower limit of the metal conduction band by extending the constant bulk potential to 
a small distance outside Zim- Our choice for Zim is a little larger than corresponding values 



obtained by fitting LDA calculations It is, however, sufficiently realistic within the 

overall precision of our model. 

The capture of electrons from a solid leads to a linear surface charge distribution along 
the surface-projected path of the projectile (Fig. |^). For a metal surface, these local charge 
densities vanish instantaneously and do not influence the charge exchange sequence or the 
motion of the projectile due to high surface plasmon frequencies of the order of lO^^sec"^. 
In contrast to metal surfaces, typical decay times for excess surface charges on insulators are 
by far too long to compensate local charge accumulations at the collision time scale. As a 
consequence the ion is followed by a linearly stretched trail of surface charges qi{xi) generated 
at times tj, i = 0, . . . , imax{Xp) at locations Xi, for which the projectile is at distances R{ti) 
above the surface. The charge depends on the ion-surface distance R through the implicit 
dependence of Xi on R. 

We assume the excess surface charges to decay exponentially with a time constant r. We 
can approximate r by having charge currents j restore the electric neutrality as depicted in 
Figure ||. The driving force of these currents are electric fields E = aj caused by the local 



surface charges 



qi{xi, t) = qi{xi) exp ( -- — - ) , t > U (2.9) 



due to local charges qi{xi) generated at times ti, i = 0, 1, 2, . . . when charge transfer took 
place. By choosing a hemisphere around the excess surface charge as Gaussian surface, 
Gauss' theorem leads to the decay time r = (27rcr)^^ in terms of the macroscopic conductivity 
cr. Typical values for a are 10^(f2cm)~^ for metals and cr ~ 1 ■ 10~^(ficm)~^ for ionic crystals, 
such as a LiF crystal, at room temperature. For LiF this simple estimate yields r ~ lO^^s 
which lies about seven orders of magnitude above the collision time of typically lO^^^s. This 
means that after capture sets in and while the HCI continues to interact with the surface, a 
positive linear charge distribution 

'■max 

X{x,t) = qi{xi)exp[-2na(t -ti)]6{x - Xi) 

i 

~ qi{xi)5{x - x.i) = \{x) (2.10) 

i 

remains on the surface-projected projectile path on the ionic crystal's surface. This charge 
distribution pulls down the potential barrier and tends to repel the HCI from the surface. 
Its contribution to the total effective potential (pTJ) amounts to 



e(0)-l \[q,,ceii] , /-^-W^^, A(x',t) 



V^oca/(x,^,t) = -^)^(-— dx' ' ^ = . (2.11) 

e(0) + 1 Z Jxo . ^2 + _ 



The term [qj^ceii] represents the integer charge withdrawn from the active surface cell by an 
active electron that crosses the barrier, e.g. [qj,ceii] = 1 for < qj^ceii < 1- The integral in 
( p.ll| ) constitutes an average over previously transferred (non-integer) charges. For metal 
surfaces, both A and Viocai vanish. 

The positive excess charge X{x')dx' within a small interval dx' near x' polarizes the 
surrounding ionic crystal. This polarization effectively screens the local charge X{x')dx' 
and is the origin of the Mott-Littleton correction |^^^ to the ground-state energy of 
the crystal. This screening correction is approximately included in ( |2.1lD in terms of the 
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static dielectric screening function (e(0) — l)/(e(0) + 1). A Mott-Littleton correction is also 
included in the work of Borisov et al. and Hagg et al. [Q. Instead of our multiplicative 
screening factor in ( |2.11| ) Borisov et al. introduce this correction as an additive contribution 
to the attractive interaction of the active electron with the left-behind hole on the surface. 
The approach of Hagg et al. resembles our approximation in that the screening of the 
capture-induced surface charge is included as a multiplicative factor which asymptotically, 
for 2 oo, becomes equal to a frequency independent dielectric screening function; it differs 
from (|2.11|) due to the inclusion of (i) dynamical screening of excess surface charges and (ii) 
fractional ionicity effects close to the surface in reference ^ . 



We note that in our version of the COM a continuous classical charge current is used to 
represent charge transfer. In our discussion of local surface charges, the discretized charge qi 
is used for convenience only and corresponds to the (small) portion of an elementary charge 
that is transferred during one timestep of the numerical propagation (see Section [II C4j , 



below). In contrast, Borisov et al. |^ and Hagg et al. [^] enforce charge quantization and 
consider local excess surface charges of at least one positive elementary charge. 



B. Local workfunction 

We now examine the dynamic change of the local workfunction while the projectile draws 
a certain amount of charge from a specific surface atom on an insulating surface. In order 
to estimate the local workfunction of an ionic crystal, we assume that target electrons are 
captured from an anion on the surface lattice. The energy necessary to remove a loosely 
bound valence electron from a surface anion can be approximated by the affinity E^^^^ of the 
free anion (3.4 eV for free F~ ions) and by adding the interactions of the detaching electron 
with all other target ions as a Madelung-background potential VMad,bg (Fig. |). This leads 
to the workfunction 

^(l^lnion) = Kind + '^Mad,bg{rlnion) ■ (2.12) 
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We evaluated VMad,bg at the mean radius r^„jo„ of the outermost ionic n£-shell {n = 2/ = 1 
for LiF), directly above the anion's lattice site, 

VMad,bg{rlnion) = 7^ ^ TT (2-13) 

where is a unit vector along the positive z axis. Since the contribution of the active surface 
anion is included via its binding energy, the active anion is exempted in the Madelung 
sum over all lattice sites Rj in (|2.13|) . Without taking screening effects into account, we 



assume Qj = —1 for all anionic charges and Qj = 1 for all cathionic charges in ( p. 131) . 



Equation (|2.12| ) yields values of 11.65 eV for LiF and 8.21 eV for KI and thus reproduces the 
experimentally determined workfunctions (12 eV for LiF and 8.2 eV for KI |^) sufficiently 
well within the overall accuracy of our approach. 

In order to include the effect of the net capture-induced local charge qj^ceu > residing 
within the active surface lattice cell j on the local workfunction, we added qj^ceii to the 
original anionic charge Qj = —1 and obtain the new charge q = Qj + qj,ceih & corresponding 
new ionic binding energy and, by using the Cowan code [^, an adjusted ionic radius 
''^ anion- Inserting these quantities into both terms of ( |2.12| ) supphes the adjusted workfunction 
iy(r^„jon)- Since the dynamical COM simulates continuous charge currents, we interpolate 
between discrete r^^j^^-values. In this way we can compute the local workfunction (p.l2|) as 
a function of the capture-induced local surface charge Qj^ceii within the active anion's unit 



cell. In contrast to the approach of Hagg et al. I^g, we do not include fractional ionicity and 
screening effects in the local workfunction. In agreement with the papers of Borisov et al. 
| 35| and Hagg et al. I^Gf, our local workfunction includes a Madelung sum for the interaction 



of the detaching/ionizing negative charge with the ionic crystal and an additional term for 
the interaction with the excess surface charge on the active anion site. However, in contrast 
to these authors, we relate the additional term to the affinity of a free anion of charge 
Qj = —1 and interpolate (using atomic ionization potentials) to effective charges g > — 1, 
as dictated by the non-charge-quantized version of the COM. For LiF, we allow for at most 
one electron to be captured from an active F~ site, such that qj^ceii < 1- For KI, we take the 
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large number of outer shell electrons on I~ in account by removing this restriction. 

In Figure | we show the various contributions to the total electronic potential for 
50 keV Xe^^^ ions that approach a LiF surface at a grazing angle of 1°. The potentials are 
displayed along an axis perpendicular to the surface that includes the projectile nucleus. The 
projectile is on the incident part of the trajectory at a distance R = 10 a.u. in front of the 
surface. The remaining projectile charge qp{R = 10) amounts to 2.8, i.e. J^*" X{x)dx = 12.2. 

C. Transition rates 

In this section, we summarize the approximations that lead to simple analytical expres- 
sions for resonant and Auger transition rates, closely following references 0,0. These rates 
are then combined in form of a system of coupled classical rate equations in order to describe 
the occupation dynamics of projectile levels. 

1. Resonant gain 

We represent the electronic structure of the projectile by its spectrum of energy lev- 
els en{R) and their occupations an{R). Both quantities change during the motion of the 
projectile. We assume hydrogenic shells with binding energies 

that depend on the effective charges 

Cleff,n{R) = (lnuc,p — X! ^ n ,n' 0"n' {R) ■ (2.15) 
n' 

The matrix Sn,n' accounts for screening effects and is determined under the simplifying 
assumption of full inner screening and no screening by outer and equivalent electrons {Sn,n' = 
1 for n > n' and Sn,n' = otherwise). The classical model of a continuous charge flow over 
the potential barrier in conjunction with discrete energy levels ( p.l4| ) requires the definition 
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of energy bins. We designate energy bins by [e„(-R)]. Each bin is attributed to a projectile 
shell n and includes electronic energies e within the interval 

(2.16) 

including corrections for level shifts due to image charge interactions and localized surface 
charges. Thus the charge current / transferred classically from the surface into the energy 
bin [£„(/?)] is considered to feed the n'th energetically shifted shell of the projectile. 

The resonant gain current (i.e. the resonant gain rate) /^^ from the surface into a par- 
ticular n-manifold of the HCI is given by the product 

rjiR,t)=a{R)UR,t) (2.17) 

of the current density 

I rmin{-W,en+i/2) 



Jn{z) = - / dED{E)p{E - H), (2.18) 

4 Jmax{V},,£n-\/2) 

and the cross section a. The conduction band density of states in free-electron-gas approxi- 
mation is given by D{E) = 'n'^\JE — Vq, where — Vq is the lower valence band limit with 
respect to the ionization threshold (cf. Fig. m). Values for the targets investigated in this 



work are Vq = 10.6, 10.9, 15.9, 14 and 10.5 eV, respectively for Au(polycrystalline) |T0|,^ 



Au(llO) Al [gSf, LiF and KI F is a volume which we assume to be one 



(a.u.)'^ for the following. The factor \j'^{E — Vj,) is the classical velocity of active electrons 
while passing the potential barrier. The energetic bottom of the conduction band lies Vq 



below the ionization threshold, and the geometrical factor 1/4 in (|2.18|) relates the isotropic 
density of states, D{E), to the electron current along the positive surface normal. The cross 
section 



a = n 



(2.19) 



2 

is equal to the classically allowed area over the potential saddle through which the current 
representing active electrons needs to flow. The effective width Ax{t) = \xi(t) — X2{t)\ of 
the saddle at any time is given implicitly by the two solutions Xi and X2 of 
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en{R) + Vim,p{qp, X - Xp, Zb + R) + Vim,e{Zb) + Vlocal{x, Zb, t) = VtoMp, X, Zb, Xp, R, t) (2.20) 



where, as in ( p.4|) , the nuclear image potential is evaluated on the axis x — Xp. 



We note that the valence band density of states of an ionic crystal per se is poorly 



represented with the free electron gas model. Hagg et al. ||36| used classical Monte Carlo 
techniques in order to simulate the over-barrier dynamics of target electrons that are released 
from an anionic center by the highly charged projectile. Their Monte Carlo study indicates 
that the electron is effectively captured from LiF at a projectile surface distance that is 
about 3 a.u. closer to the surface than the onset of classical electron capture across the 
potential barrier. Interestingly, our numerical results, including capture rates modeled with 
reference to the free electron density of states, also indicate a delayed onset between the 
initiation of a classical over-barrier current and the projectile-surface distance where one 
elementary charge has been transferred (see Section IllA below). Furthermore, we point 
out that previous dynamical COM studies on electron capture in collisions with Ceo [0 have 
shown a rather weak dependence of final projectile charge states and critical capture radii 
on variations in the resonant capture rates. We therefore conclude, that within the overall 
accuracy of the COM and in view of the narrow valence band of ionic crystals, the free 
electron gas model is sufficiently realistic for providing acceptable estimates for resonant 
capture rates. The agreement of our simulated projectile energy gains with experiments for 
ionic crystal surfaces provides further support for this approximation (see Section IIIB2, 
below) . 

2. Resonant loss 

The rate of electron loss from atomic energy levels into unoccupied bulk levels can be 
obtained from the electron's orbital frequency of revolution 

/n = 1^ (2.21) 

multiplied with its probability to hit the saddle region that is approximately given by 
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P(e„) = [ZcritiSn) - Z]/Z„it{£n) (2.22) 

where Zcriti^n) denotes the critical distance where the first electron capture occurs into shell 
n i, 

i:^{z) = U-P{e^). (2.23) 

For insulators, particularly for LiF with no unoccupied band levels below the vacuum level, 
resonant loss is irrelevant, and Jj^' = 0. 

3. Auger processes 

Intraatomic Auger transitions induce small changes in the projectile occupation during 
the interaction with the surface. Following reference |Q we express the Auger rates by a 
simple analytic fit through data points 



5.06-10-3 



that have been calculated with the Cowan code |57] for fast transitions between two given 
shells rii and nj. 

4- Rate equations 

The dynamically varying projectile populations are obtained as solutions to the coupled 
set of classical rate equations 

^ = Q{An - an)Tlf - a„r^' + ^ r„/,„Wj,„/ - 2Wi^n ^n,n'Wfy (2.25) 

n'>n n'<n 

where the degeneracy of shell n is given by An = 2n^. is the unit step function. The 
(empirical) statistical factor Wf^n = l/(l + 1.5an) corrects for the decrease in Auger transition 
rates due to increasing populations a„ of the final level. The statistical factor Wi^n = |an(ctn — 
1) takes the equivalence of electrons in the initial shell into account. 
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D. Projectile motion 



Before reaching the first critical over-barrier radius, the motion of the incoming HCI is 
solely affected by its attractive self-image force. For metal surfaces this force is given by 

(^' = -(2(7^)' P-^'" 

and for insulating surfaces by the derivative of ( |2.4| ). After the projectile has reached the 
first critical over-barrier distance, charge transfer begins and, for insulators, the self-image 
force ( |2.26| ) starts to compete with the repulsive force created by localized surface charges 
( p.ll| ). For grazing collisions, the latter force is weak, due to the large projectile velocity 
component parallel to the surface, which rapidly increases the distance between previously 
created surface charges and the HCI. 

At distances R smaller than the largest radius of occupied atomic orbitals, (r)„, the 
electron clouds of the incoming HCI and the surface ions begin to penetrate each other. The 
accurate description of this situation would require detailed quantum dynamical calculations, 
which are far beyond the overall simplistic nature and accuracy of the COM. In order to 
determine the classical motion of the projectile, we employ the Thomas-Fermi model and 
use the Thomas-Fermi-Moliere potential energy 



TI/ / s <lnUC,p ■ (lnUC,t 1 \ /o r)'7^ 

WTFM{r) = - • (2.27) 



This interatomic potential includes the Coulomb repulsion between the two nuclear charges 
(lnuc,p and qnuc,t of the HCI and a target atom, respectively, and a screening function (p that 
depends on the internuclear distance r scaled by the screening length a, 

=Eaiexpf-A-) (2.28) 
a ~^ \ a/ 



a = 0.88534/ ■y' qnuc,p + qnuc,t , 
with {a,} = {0.35, 0.55, 0.10}, {A} = {0.3,1.2,6.0}. 
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For for small perpendicular projectile velocity components, as typically given in grazing 
incidence collisions, the interaction of the HCI with an array of surface atoms or ions is well 
represented by taking the planar average of ( |2.27| ) |5T[] , 

^Qnuc,pQnuc,t il f ^ 



Wtfm{R) 



' -)=E^«^pf-A-) (2.29) 



a ■ '~[l3i V a J 

where, for simphcity, we have assumed a square lattice with lattice constant d. 

For ionic crystals that are composed of two different ion species, we apply ( p.29| ) sepa- 
rately to surface lattices of anions and cathions. This results in the planar averaged potential 

Wtfm{R) = W^tuiR) + "(^) (2.30) 



where W^^^f + W^'^p'l}"" are constructed according to ( |2.29| ) with g„„c,t replaced by the 



respective nuclear charges of anions and cathions and with the distance d between anions 
or between cathions, respectively. 

The force exerted on the projectile by the capture-induced surface charge distribution 
( p.lOp is repulsive with a parallel component that accelerates the projectile in positive x- 



direction. It is given by 

F (X m-'M^n f^-,,. A«t) {R,X,-x') 

The factor (e(0) — 1)/ (e(0) + 1) accounts for the static dielectric screening (cf. equ. (|2.11|) ). 
The effect of this force on the projectile motion parallel to the surface is small. For 30 keV 
;^gi5+ incident under a grazing angle of 1° on LiF, it changes the parallel component of the 
projectile velocity by 0.2 %. 

The net force on the projectile is now given by the negative gradient of Wtfm, the 
projectile self-image force (p.26|) , and ( p.31| ). 



Fnet{Xp, R) = I — — (WTFAf(-R)) + Fim,p-p{R) | 6^ + Fiocal{Xp, R). (2.32) 
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III. NUMERICAL RESULTS AND DISCUSSION 



In this section we are going to discuss the pronounced differences in the above-surface 
neutrahzation dynamics between insulating and conducting targets as predicted by our sim- 
ulations. All electronic interactions depend strongly on the projectile's position on its tra- 
jectory {Xp{t), R(t)), which, in turn, depends on the charge state evolution qp{t) and, for 
insulator targets, on the trail of positive excess surface charges. Due to this coupling of 
nuclear and electronic degrees of freedom, the reconstruction of measured projectile deflec- 
tion angles and projectile energy gains will supply support for the specific interaction model 
implemented in our simulation. Even though it is not yet possible to extract direct evidence 
from the experimental data for the time evolution of many quantities included in our simu- 



lation, such as shell populations and level shifts (Sec. [Ill A|) , good agreement between energy 



gain measurements and theory (Sec. IIIB) supports the validity of our model assumptions. 



A. Interaction dynamics 

Fig. ^ shows the simulated projectile charge evolution for Xe^^^ ions colliding with an 
Al target in one case, and a LiF crystal in the other. The incident projectile energy is 
50 keV at a grazing incidence angle of 1°. For Al, the first critical distance for classical 
over-barrier capture is Rc ~ 38. In the case of LiF, the interplay of the large workfunction 
and image charges, which compared with the Al target are reduced by the altered dielectric 
response function, effectively shifts the onset of charge transfer by about 10 a.u. to Rc ~ 30. 
Our version of the COM does not impose charge quantization. By rounding to nearest 
integer charges, the neutralization sequence is therefore completed when the projectile charge 
becomes smaller than 0.5 . In comparison with LiF, we find that for the same ion-surface 
distance R the early onset of electron transfer on Al leads to smaller projectile charges above 
the metal target. 

Differences in the time evolution of resonant gain processes become apparent by com- 
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paring resonant gain rates (Fig. Pd and Fig. |^) or projectile level occupations (Fig. |c and 
Fig. for the two targets and incident Xe^^"*" ions. On the Al target, the n=18 shell of 
the projectile becomes strongly populated with filling rates of the order of lO^^s"^ (Fig. 0). 
For LiF, resonant gain transfer occurs into the n=14, 13, and 12 shells. The sharp maxima 
of the gain rates on LiF exceed the metal rates by almost two orders of magnitude. The 
average neutralization rates, however, i.e. the slopes of the corresponding shell occupations 
(Fig. Pa and (Fig. 0b), are only slightly higher in LiF. The interruption of the projectile 
neutralization between R = 19 and R = 21 coincides with the transient increase of Vb above 
W in Fig. ^(b). Resonant loss processes are either forbidden (LiF) or contribute with neg- 
ligible rates (Al). For the Al target, shells below the resonantly populated level n=18 are 
populated in Auger transitions. 

The regularly spaced spikes in the resonant gain rate for Xe^^"*" impinging on LiF (Fig. ^) 
originate in the capture-induced local surface charges. As the surface-projected path of the 
ion enters a new surface cell containing a single fluorine F~ ion, the workfunction is reset 
to its original value W = 12 eV, such that the local Fermi level suddenly moves upwards 
thereby stimulating over-barrier capture. The corresponding workfunction changes amount 
to up to 6 eV (Fig. |^), whereas oscillations in the barrier height H of the potential saddle 
remain comparatively small with an amplitude of less than 1 eV as the ion travels over the 
surface cell boundaries. This relatively inert behavior of VJ, can be explained by the moderate 
influence of the local surface charges on the total potential Vtot near the saddle position Z}, 
which is situated typically a few a.u. in front of the first bulk layer (see also Fig. We 
note that, if capture from a given anion proceeded, the local workfunction would increase 
by about 10 eV per unit capture-induced surface charge. However, due to the high transfer 
rates of up to 10^^s~^ the continuous current of negative charge is quickly cut off at the 
moment when the Fermi level is shifted below the saddle point Vf, ("over-barrier cut-off"). 
In other words, the local workfunction change (i.e. the shift of the Fermi level towards lower 
energies and below Vf,) generated by capture-induced surface charges produces the peaked 
structures in the resonant gain rates. For the Al target, the absence of local surface charges 
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results in a comparatively steady evolution of the gap between Vb and W (Fig. ^ja) and 
results in the mostly smooth development of the dominant resonant gain rates in Figure |^a. 

Considering the characteristic discrepancies in the resonant exchange mechanisms, it is 
surprising that the average rate of neutralization is very similar for both targets. The effects 
of the low alkali halide workfunction on the onset of charge exchange, the reduced dielectric 
response of the insulator, and local surface charges appear to be counterproductive. As will 
be shown below (Sec. [1IIB|) , this interplay is also related to the strikingly small differences 



between the image energy gains of a particular HCI on LiF and Al targets. 

In Figure ^ we compare ion trajectories for grazingly incident 20 keV Xe^^^ on Al and 
LiF surfaces. At large distances, the magnitude of the perpendicular velocity component \vz\ 
steadily increases due to the attractive projectile self-image force. The short range TFM- 



potential in ( p.27|) causes the inversion of the trajectory in a small region that measures 
about 2 a.u. relative to the vertex of the trajectory leading to nearly specular reflection. 
For the LiF target the attenuated dielectric response of the insulator weakens the image 
attraction in comparison with a metal target. We simulated this effect in a separate calcu- 
lation for LiF where we replaced the insulator specific dielectric response in the projectile 



(self-) image interactions by the asymptotic response of a perfect metal, taking ( p^.26D and 



the limit e ^ oo in (|2.4|) . This yields a noticeable change in the ion trajectory. Before 
capture sets in at large R, only the image force acts on the HCI and the "metal dielectric 
response" trajectory for LiF nearly coincides with the ion trajectory in front of Al. The re- 
placement of the insulator-specific dielectric response by the metallic response on LiF moves 
the onset of charge transfer 13 a.u. closer to the surface. This shifts the potential barrier 
upwards due to the more repulsive (unscreened) projectile image term ( p.4|) , thus reduces 
the amount of charge captured, and more than doubles the overall energy gain above the 
target, clearly leading away from both our dynamical COM with insulator-specific response 
and experimental results (see Sec. [ill B| ). 

In another separate simulation we have eliminated all effects due to the capture-induced 
positive surface charge distribution on LiF (Fig. The local surface charges add to the 
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projectile repulsion near the surface. The vertex is now located about 0.1 a.u. closer to the 
first bulk layer than for the genuine 50 keV Xe^^"*"-LiF simulation including local surface 
charges. This small shift suggests that the direct influence of these surface charges on the 
projectile trajectory is rather small. The large discrepancy of 12.8 eV in the image energy 
gains between both simulations (see Section [111 B 21 ) originates mainly from deviations in 
the neutralization dynamics leading to a higher average projectile charge in front of the 
surface when local surface charges are disabled. This can be understood by considering 
that the component Viocai in (|2.11|) pulls down the potential barrier ( p.l|) and also the 
projectile energy levels ( p.l6|) and thus counteracts electron loss due to level promotion 
into the continuum (if an occupied level gets promoted to the continuum, we assume that 
the level is instantaneously ionized). 

The distance of closest approach to the surface under grazing incidence is determined by 
the initial velocity component perpendicular to the surface, the total kinetic energy gain of 
the HCI at the turning point, and by the composition of the target material via the screened 
interatomic interaction ( ^.271 ). Our simulations yield turning points at -R = 1.1 and 1.3 for 
Xe^^^ ions impinging on Al and LiF, respectively, at an incident energy of 30 keV and an 
incidence grazing angle of 1°. 



B. Image energy gains 

After having presented detailed results on the interaction dynamics in the previous sec- 
tion, we are now going to demonstrate that the extended dynamical COM can quite ac- 
curately reproduce previously published (measured and simulated) data on image energy 
gains for both conducting and insulating crystals over a wide range of initial projectile 



charge states. Our simulations as well as recent experiments show that the neutraliza- 
tion of the HCI is completed prior its reflection for a wide range of initial projectile charge 
states. The inversion of the perpendicular velocity component v± takes place at a distance 
of a few atomic units above the surface (Fig. ^, and the measurable difference between the 
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asymptotic incident and reflection angles of an ion beam can be straightforwardly correlated 
to the net image energy gain ||5^j27[] . 



1. Metals 

Image energy gains of a HCI impinging on metal surfaces are characterized by an ap- 
proximate gp/^-increase with the initial projectile charge state P, p!0| , p7| . A lower limit 
for the energy gain can be deduced by assuming that the projectile is instantaneously and 
completely neutralized at the first critical distance Rc ~ ^f^Qp + ^/ (2W^) . We shall 
refer to this estimate as "simple COM". The energy gain for large Qp is then given by the 
analytical formula 

Consequently, /W should be independent of the target material. 

The simple COM can be improved by letting one electron transfer to the HCI each time 
the over-barrier condition is fulfilled at consecutive critical radii for the first, second, etc. 



capture. This version of the COM is called the "staircase model" |T^. In contrast to the 



dynamical COM (Sec. [11 A] ), in the simple and staircase model the charge transfer current is 



quantized. For energy gains of Xe''+ projectiles on an Al surface, the staircase model almost 
coincides for all initial charges Qp with the dynamical COM, and the simple model predicts, as 
expected, lower energy gains for all qp (Fig. 0). Except for the highest charge states, both, 
the staircase and dynamical COMs agree with the experimental gains of Winter et al. pl, 



even though the more elaborate dynamical COM employs transition rates that depend on 
the width and depth of the potential saddle. The simple model underestimates the measured 
energy gain, except for the highest charge states, where agreement with the experiment may 
be fortuitous. Except for the simple COM, all simulations intersect the experimental error 
bars for charge states < 30. All simulations show the general gp^^-trend. 

The deviation in the experimental data from the approximate gp/^-proportionality of the 
energy gain in all COM versions above qp ~ 26 (Fig. |l^) has been scrutinized by Lemell et al. 
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10| . The authors rule out both saturation effects in the surface charge density fluctuations 
induced by the HCI a.t R = Rc and effects due to the parallel velocity of the HCI. They 
conclude that the measured deviation from the gp/^-proportionality is due to incomplete 
screening of outer shells at decreasing R which, for high initial projectile charges, leads to a 
faster decrease of the effective projectile charge and, therefore, to a diminished increase in 
the energy gain as a function of qp. However, as far as we know, the initial charge state qp 
at which the experimentally observed plateau appears has not yet been reliably calculated 



within any COM (see also the review article of Winter and references therein). Figure ^ 



also shows that the staircase COM calculation of Lemell et al. ||T0[ agrees with our results. 

Kurz et al. have analyzed total electron yields for higher charge states as a function of 
the inverse projectile velocity. Their data for Xe''"'', g = 34 ... 50 and Th'^'*', g = 61 ... 79 on 
gold surfaces under perpendicular incidence provide, if at all, weak evidence for a deviation 
of the energy gain from the g^/^-proportionality (Fig. |ll]). We note that the experimental 



method of reference |^4[ is prone to larger errors than the deflection angle method |Q . Our 
dynamical COM data for Thorium are near the upper end of the experimental error bars. 

In Figure we compare our dynamical COM results for 150 keV ions on Au with energy 
gains measured by Meyer et al. [^] and with the COM simulation of Lemell et al. []IU|. The 



experimental results show good overall agreement with our calculations for both projectiles 
but fall systematically short of the dynamical COM values above gp ^ 30. 

2. Insulators 

Energy gains for 50 keV Xe ions directed under a grazing incidence angle of 1° on alkali 



halide crystals (LiF and KI) have been measured recently by Auth et al. []3T| , |33| by using 



the deflection angle method ||52|. Our extended dynamical COM simulations agree with 
experiment for the Kl target for qp < 17 (Fig. |13|). However, experiment and simulation 
tend to deviate in a systematic way for the LiF target, where, for low and intermediate 
incident charge states, the measured values slightly exceed our simulations. We tried to 
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identify an adjustable parameter in order to further improve the agreement with experiment 
for both targets. In the oncoming paragraphs we will discuss several of the effects that 
appear in our insulator extension of the COM and use Xe^^"*" on LiF and KI, with an energy 
gains of 43.6 eV and 31.8 eV, clS 8b reference. 

At first we take a closer look at effects that are induced by the surface charges. We 
observe that the restriction of one removable charge per F~ ion (which we did not apply 
to iodine for its vast number of outer shell electrons) lowers the energy gain by 1.6 eV. 
Furthermore, disregarding all capture-induced surface charges (cf. curves labeled qiocai = 
in Fig. |l^) increases the energy gain of Xe^^^ on LiF to 56.3 eV. For the KI target, however, 
the neglect of surface charges increases the energy gain to 37.0 eV for incident Xe^^^ which 
lies above the experimental error bars. Our numerical results show the expected increase 
for energy gains at all charge states if we discard surface charges. For the KI target, the 
inclusion of capture-induced surface charges improves the agreement between simulated and 
measured energy gains. 

With respect to the ionic conductivities, we note that a has to be increased by more than 
six orders of magnitude in order to induce any significant change in the energy shifts. Despite 
our crude estimates for a and the order-of-magnitude derivation of the time constant r in 
( ^^ , we can therefore exclude life-time effects of capture-induced surface charges on the 
simulated energy gains. 

The image plane is located at one anionic radius above the uppermost bulk layer. This 
choice constitutes an upper limit for Zim- Alternatively, as an accurate value is difficult to 
assess and the concept of an image plane is not well defined for ionic crystals, one could use 

as an adjustable parameter. Placing the charge distribution at the topmost lattice plane 
diminishes the image energy for Xe^^^ by 1.6 eV on LiF. For KI the energy gain slightly 
increases by 0.6 eV. 

We have performed all simulations with free-electron densities of states and a constant 
volume factor F = 1 in the transition rates ( |2.17| ), which is a rather poor approximation for 



an ionic crystal and may not give sufficient credit to the characteristics of a particular crystal; 
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the negative fluorine ion possesses six 2p-electrons, whereas iodine ion holds a large number 
of loosely bound electrons. In an attempt to work this information into the simulation, we 
reduced the resonant gain rates ( p.lT] ) for the LiF target by a factor of 5. As a result we find 
that these modified rates lead to energy gains on LiF that lie inside the experimental error 
bars for all charge states. In a more realistic representation of the target electronic structure, 
more attention must be given to the valence electrons of the anions. The above-mentioned 
discrepancies for KI and higher charges of the incident projectile may be related to the 
simplified representation of the target electronic structure inherent in our implementation 
of the dynamic COM. 

As explained in Section IIA, our simulations were limited to trajectories with 1^ = 0, 
for which the collision plane intersects anionic and cathionic nuclei along the [100] direction. 
Since in surface scattering experiments the incident ion beam illuminates a surface area 
that is large compared with a surface unit cell, we addressed the sensitivity of our simulated 
energy gains to changes in Yp. For Yp = 0.5d, corresponding to surface-projected trajectories 
half way between ionic rows, we find for 30 keV Xe^^"*" projectiles incident under 1° on LiF 
a kinetic energy gain of 44.6 eV, compared to 43.6 eV for trajectories with Yp = and 
an experimental value of 53.2 eV. The slightly larger energy gain is consistent with 
the increased average distance of capture-induced surface charges from the projectile. The 
change in energy gain with Yp is sufficiently small such that, within the overall accuracy of 
our calculation, we do not need to include a time-consuming average over trajectories with 
different Yp in our simulation. 

IV. CONCLUSIONS 

In this work we applied and discussed extensions to the classical over-barrier model 
that include insulator specific effects such as capture-induced local surface charges, local 
workfunction changes, and the dielectric response of the target. A detailed study of the 
interaction mechanism has been presented in terms of the time evolution of projectile level 
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occupations, transition rates, and several other quantities involved in the neutralization 
process. 

Our results are in good agreement with previously published experimental data for highly 
charged ions impinging on two different alkali halide ionic crystals. In order to verify the 
basic framework of our implementation of the dynamical COM, we have disabled all effects 
related to insulators and found good agreement with energy gain measurements for a variety 
of incident ion charge states and metal targets. 
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FIGURES 

FIG. 1. The band structure of the ionic LiF crystal with a large workfunction of 12 eV and 

a wide band gap of 14 cV. The (polycrystalline) gold target represents a typical metal with a 
workfunction of about 5 eV and a continuum of unoccupied conduction band states above the 
Fermi level. 

FIG. 2. The linear capture-induced charge distribution X{x,t) trailing the path of the ion 
(schematically). The capture sequence starts when the incident highly charged ion reaches the 
critical over-barrier distance at a position (Xp,R) = (xo,Rc) at time tQ. For our applications to 
ionic crystals, we assume that electrons are captured from the closest surface anion. 

FIG. 3. Classical currents j that are driven by the field E = aj of a capture-induced surface 
charge restore electric neutrality. A decay time constant r can be derived from the macroscopic 
conductivity a by applying Gauss' theorem to the current j and a Gaussian surface given by a 
hemisphere with the positive excess surface charge in its center. 

FIG. 4. The charge state dependent local workfunction W of LiF is approximated by splitting 
potentials acting on a bulk atom into an anionic contribution, the ionic binding energy E^-,^^^, and 
the Madelung background potential VMad,bg representing the rest of the bulk. The distance r^^nicm 
denotes the orbital radius of the most loosely bound subshell of the ionic state and is dynamically 
adjusted to the excess local capture-induced charge q. 

FIG. 5. Contributions to the total electronic potential Vtot for 50 keV Xe^^"*" approaching a 
LiF surface at an angle of 1° along an axis perpendicular to the surface that includes the projectile 
nucleus. The projectile is on the incident part of the trajectory at i? = 10 and qp{R = 10) = 2.8. 

FIG. 6. Results for Xe^^"*" at E^in = 50 keV and an incidence angle of 1°. Time evolution of 
the projectile charge state for LiF and Al surfaces (a). Resonant gain rates on LiF for the two 
highest resonantly populated shells (b). Projectile shell occupations on LiF (c). 



30 



FIG. 7. Results for Xe^^+ {Ekin = 50 keV) ions impinging under 1° grazing incidence conditions 
on an Al surface. Time evolution of the resonant gain rates for the highest resonantly populated 
shells (a). Shell occupation (b). Charge exchange primarily takes place via resonant gain into the 
n = 18 shell. 

FIG. 8. Results for Xe^^^ ions impinging with i?fcm = 50 keV at an grazing angle of 1° on 
Al and LiF surfaces. The plot shows the time evolution of the potential barrier T4, the target 
workfunction W, and the projectile shell occupations of the most active shells. Vf, and W display 
characteristic oscillations on LiF, whereas all potentials evolve smoothly on Al (see text for details). 

FIG. 9. Projectile trajectory in terms of the perpendicular projectile velocity component 
versus R for 20 keV-Xe^^"^ ions impinging under 1° on Al and LiF. The first two curves exhibit 
standard dynamical COM simulations on these targets. The next two curves show simulation 
results for the same Xe^^"*" projectile on LiF when local surface charges have been disabled in the 
third and a large "metal" value for the dielectric susceptibility e oo has been chosen for the 
fourth curve. Xp = corresponds to the vertex of the full simulation for LiF. 



FIG. 10. Experimental [^], simulated staircase COM results ||l^, and our simulated energy 
gains (dynamical COM) for Xe'^^ (3.7 q keV, 1.5°) on an Al surface. The simple model assumes 
instantaneous complete neutralization at the first critical distance Rc and sets a lower boundary 
for projectile energy gains. The staircase COM instantaneously transfers one charge unit each 
time the over-barrier condition is fulfilled. In the dynamical COM continuous charge currents flow 
between projectile and surface with rates derived from a classical model. 



FIG. 11. Experimental and our simulated data (dynamical COM) for very high charge 
state ions impinging on polycrystalline gold. 



FIG. 12. Experimental [28|, simulated staircase COM results |1C], and our simulated data 



(dynamical COM) for 150 keV and Pb'^^ ions with charge states qp < 36 on Au. 
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FIG. 13. Experimental energy gains |^3| compared with our dynamical COM simulations, 
suits obtained by neglecting capture-induced local surface charges are labeled as qiocai = 0. 



Table 1. Static limit (eo), optical limit (eoo) and characteristic frequency ujq used 
for the dielectric response of LiF and Kl crystals at T=290°C 





eo 


Coo 


ujo in 10 ^ a.u. 


LiF 


9.00 


1.93 


1.39 


KI 


5.09 


2.65 


0.46 
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